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ABSTRACT 

Infrared dark clouds (IRDCs) are dense, molecular structures in the interstellar medium that can harbour sites of high-mass star 
formation. IRDCs contain supersonic turbulence, which is expected to generate shocks that locally heat pockets of gas within the 
clouds. We present observations of the CO J = 8-7, 9-8, and 10-9 transitions, taken with the Herschel Space Observatory, towards 
four dense, starless clumps within IRDCs (Cl in G028.37+00.07, FI and F2 in G034.43+0007, and G2 in G034.77-0.55). We detect 
the CO J = 8-7 and 9-8 transitions towards three of the clumps (Cl, FI, and F2) at intensity levels greater than expected from 
photodissociation region (PDR) models. The average ratio of the 8-7 to 9-8 lines is also found to be between 1.6 and 2.6 in the three 
clumps with detections, significantly smaller than expected from PDR models. These low line ratios and large line intensities strongly 
suggest that the Cl, FI, and F2 clumps contain a hot gas component not accounted for by standard PDR models. Such a hot gas 
component could be generated by turbulence dissipating in low velocity shocks. 
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1. INTRODUCTION 


Giant molecular clouds contain significant supersonic turbu¬ 
lent motions and this turbulence may play a key role in sup¬ 
porting molecular clouds against collapse on large scales, as 
well as in generating small-scale density perturbations that can 
then locally collapse to form stars (e.g., iPadoan & Nordlundl 
2002tlMcKee & Qstrikeiil2007h . Numerical simulations of mag¬ 
netohydrodynamic (MHD) turbulence, similar to that ob¬ 
served in molecular clouds, show that such turbulence de¬ 
cays on the order of the c rossing time at the driving scale 
’ Gammie_<&^9striker 1 9961: iM ac Low et al.l 119981 : IStone et aT 


1998; Mac Lovd 1999 •iPadoan & Nordlund 1999; Ostriker etal 


20011) . In typical Galactic star-forming regions, the volume av¬ 
eraged turbulent dissi patio n rate is on the order of the cosmic 
ray heating rate dPon et alJl2012 ). This turbulent heating, how¬ 
ever, only occurs within localized shock fronts, such that these 
shocks create a hot gas component with a low volume filling 
factor, thereby producing a very different spectral signature than 
would be expected from uniformly heated gas (iPon et alJl2012l) . 

iPon et al] (1201 2l) present MHD, C-type shock models for 2 
and 3 km s 1 shocks propagating into gas with densities between 
10 2 5 and 10 3 5 crcT 3 . By scaling these models to the expected 
turbulent energy dissipation rate of a molecular cloud (e.g.. 


* Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


iBasu & Muralil 120011) . IPon et al.l (1201 2l) predict the integrated 
intensities of shock excited lines from low-mass star-forming 
regions. They show that most of the energy in these shocks 
is radiated by rotational transitions of CO and, by comparing 
their shock results to the photodissociation region (PDR) mod¬ 
els of iKaufman et al.l (1 9991) . they predict that mid-/ CO lines 
(Jupper > 6) should be dominated by emission from shocked 
gas, despite this hot shocked gas only having volume filling fac¬ 
tors of the order of 0.1%. Such excess emission in the CO J 
= 6 —> 5 line, consistent with an extra hot gas component, has 
already been detected in the Perseus Bl-East 5 low-mass star¬ 
forming region jPonetal.1120141) and in the Taurus molecular 
cloud (iLarson et al.li2015 ). 

Infrared dark clouds (IRDCs) are molecular gas structures, 
which are usually filamentary, that appear as dark lanes against 
the bright mid-infrared (~ 8/im) Galactic background due to 
their high column densities, with visual extinctions_approach- 
ing 100 mag in s ome regions (e.g., iRathborne et al.l [2006; 
iButler & Tan 120121) . Some IRDCs are associated with high- 
mass star-forming regions. Four particularly dense, infrared dark 
clumps within IRDCs (Cl in G028.37+00.07, FI and F2 in 
G034.43+00.24, and G2 in G034.77-00.55) were selected from 
the lButler & Tanl (120121) sample as prime candidates for evolved 
clumps without any embedded young stellar objects. These four 
clumps are ideal candidates to search for a small volume filling 
factor, hot gas component caused by the dissipation of turbu¬ 
lence, since they are unlikely to contain significant internal heat- 
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ing sources that can mimic the heating caused by low velocity 
shocks and because they have large columns of gas in which the 
shocks can exist. 

The four clumps observed are described in more detail in 
Sect. In Sect. [3] observations taken with the Herschel Space 
Observatory of the 12 CO / = 8 —» 7, 9 —> 8, and 10 —> 9 transi¬ 
tions towards the Cl, FI, F2, and G2 quiescent clumps are pre¬ 
sented. Photodissociation region (PDR) models are compared to 
the observations in Sect. [4] The implications of such a compari¬ 
son are discussed in Sect.[5]and we summarize our main findings 
in Sect. [6] 

2. SOURCES 

We consider clumps to be structures that are large enough to 
fragment into star clusters and cores to be structures likely to 
form a single star. To match such a theoretically based classifica¬ 
tion system to observed structures in IRDCs, we choose to refer 
to objects on the size scale of a parsec as clumps, and objects 
on the size scale of 0.1 pc as cores. Such a naming convention 
helps us to clarify the size scale of the object being described 
and to differentiate between substructures within regions. With 
this working, observational definition, clumps will typically have 
masses on the order of hundreds of solar masses (IButler & Tanl 
2()09. l20T2l) while cores will have masses on the order of tens of 
solar masses dTan et al .1120 IButler et al~ll20l4 . 

The four clumps examined in this paper are the clumps stud¬ 
ied by iTan et al.l ( 201 3f). The locations of these four clumps, 
as given by IButler & Tanl (1201 2h . are listed in Table [T] and the 

clump s will be referred to as Cl, FI, F2, and G2 (Butler & Tanl . _ 

120091 120121) . These clumps are spread between three differ- Infrared (HIFI; de Graauw et al 
ent IRDCs, the C (G028.3 7+00.07), F (G034.4 3+00.24), and G 
(G034.77-0.55) IRDCs of IButler & Tanl d2009l) . IRDC C is also 
known as the Dragon Nebula ( Wandl2015l) . These clumps were 


of Rathborne et al. ( 200f j) and the 24 micron so urce towards the 
edge of the FI clump dChambers et al.l l2009h . which are dis¬ 
cussed further in Sect. [5] 

The MM7 clump, at ( 18 A :53 m 08.3 J , 1°:27’:13”) and first 
named bv lRathborne et ail (120061) , wa s previously called t he F2 
clump by Butl er & Tanl ( 2009 ). When IButler & Tan (12012) up- 
dated the central locations of 42 of the iRathborne et al.l (12006) 
clumps by selecting the highest m ass s urface density peak within 
the clump areas defined by Rathbomeetak[( 20061), the F2 clump 
was the only clump for which IButler & Tan (1201 2l) were unable 
to find a mass surface density peak within the clump bound¬ 
ary or within a few arcseconds of the boundary. The peak that 
IButler & Tan (2012|) ended u p assi gning to F2 is ~ 10 arcsec out¬ 
side of the Rath borne et al. (2006) boundary and 23 arcsec from 
the central location of MM7. Due to this position shift, we refer 
to the 0.5 pc scale mass surface density structure associated with 
the IButler & Tanl (l2012l> F2 loca tion as the F2 clump an d refer to 
the structure associated with the Rathborne et al. (2006) location 
as the MM7 clump. A closer examination of the Butler & Tanl 
( 12012 !) mass surface density map suggests that the F2 clump and 
MM7 clump are indeed different clumps (see Fig.[l}. 


3. HERSCHEL OBSERVATIONS 

3.1. Reduction 

As part of the open time project OT2_pcaselli_8, roughly 1 
arcminute square maps of the CO / = 8 —> 7, 9 —» 8, and 
10 —> 9 line emission from the Cl, FI, F2, and G2 clumps 
were obtained using the Heterodyne Instrument for the Far- 


Space Observatory (iPilbratt et alJ 


2010) on board the Herschel 


2010). All of the observations 


selected by ITan et al.l ( 20131 ) for five reasons. First, the clumps 


(Rathborne et al. 2006 

Chambers et al. 2009), mid-infrared ex- 

tinction (Butler & Tan 

2009, 

2012b, and molecular line emis- 


sion, including, but not limited to, N 2 FR and N 2 D + emission 


ng, 

dFontani et al.ll201 ll) . Second, the clump s are dark at 24 an d 70 
microns in Spitzer MIPSGAL images (ICarev et al.l I2009I) . in¬ 
dicating the lack of massive embedded protostars. Third, the 
clumps are in relatively quiescent environments, indicated by the 
lack of mid-infrared sources within 20 arcsec (Beniamin et al.l 
120031: ITan et al! 120131). Fourth, the clumps have large mass sur¬ 
face densities ( Butler & Tanl 120091 1201 2l) and fifth, they have 
large deuterium fractions measured by the N 2 D + to N 2 FR col¬ 
umn density ratio dFontani et al.l 1201 lb . These four clumps are 
thus prime candidates to be sites of star formation in the near 
future, but are unlikely to currently contain massive protostars. 

ITan et al] (12013l) detect small, 0.1 pc sized cores in N 2 D + 
emission within each of these clumps, with the Cl and G2 
clumps harbouring two local maxima of NgD + emission each. 
These individual N 2 D + cores are denoted as Cl-N, Cl-S, G2-N, 
and G2-S. 

Figure|T]shows the mass surface densities of the four clumps, 
as derived from 8 micron extinction mapping (IButler & Tanl 
2012 ), and the central locations of the six cores. The contours 
start at 0.075 g cm~ 2 (Ay of 17 mag) and increase by increments 
of 0.075 g cm 4 , with the typical uncertainty being on the or- 
der of 30% based upon unc ertainties in the opacity per unit mass 
(IButler & Tanl 120091120121) . The uncertainty base d upon back¬ 
groun d fluctuations is of the order of 0.013 g cm 2 (IButler & Tanl 
120121) . This figure also shows the locations of the MM7 clump 


were conducted in position switching mode with on-the-fly map¬ 
ping. The off positions, in J2000 right ascension and declina¬ 
tion, used for Cl and G2 were (18 ,! :42 m :28.18 f , -4°:07':22") 
and (18 / ':56"':32.14 I , 1°:24':45"), respectively, while a common 
offset position of (18 /, :53" , :26.18T 1°:26':16") was used for FI 
and F2. HIFI obtains data using both a Wide Band Spectrome¬ 
ter (WBS) and High Resolution Spectrometer (HRS). Since the 
WBS data have lower noise and the observed lines are quite 
weak, we will hereafter only use the WBS data. The noise in 
the WBS and HRS spectra are correlated, such that averaging 
the two spectra together would not produce any significant gain 
in the signal to noise of the lines. Further details of the observa¬ 
tions, including the map centres, source times, beam sizes, obser¬ 
vation bands, upper energy levels, rest frequencies, observation 
IDs, main beam efficiencies, and observation dates are given in 
Table 0 

All of the Herschel data were reduced and processed using 
the Herschel Interactive Processing Environment (HIPE). Dur¬ 
ing this reduction, multiple versions of HIPE between version 
8 and 12.1.0 were used. The Herschel data were run through 
the default pipeline before being downloaded from the Herschel 
Science Archive. Standing waves were removed from the data 
and a second order polynomial was fit to the baseline. The H 
and V polarization data were combined and spectral cubes were 
created from the data. When making the spectral cubes, the do- 
Gridding task was initially allowed to choose the optimal map 
grid for the data. After examining the data to determine which 
maps show detections, the doGridding script was then rerun, on 
a clump by clump basis, on the original data to grid all of the 
data sets showing detections onto the grid map for the highest 
J line with definitive detections. Due to the similarities in beam 
sizes for the different lines, the higher J maps were not further 


Article number, page 2 of[B] 

























































































































A. Pon et al.: Mid-./ CO shock tracing observations of infrared dark clouds I 


Table 1: Target Locations 


Clump 

Cloud 

RA (J2000) 

Dec(J2000) 

/ 

b 

Ulsr 

d 



(h:m:s) 

(°:':") 

(°) 

O 

(km s _1 ) 

(kpc) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Cl (MM9) 

G028.37+00.07 (C) 

18:42:47.0 

-04:04:09 

28.32450 

0.06655 

78.6 

5.0 

FI (MM8) 

G034.43+00.24 (F) 

18:53:16.6 

01:26:10 

34.41950 

0.24583 

57.1 

3.7 

F2 

G034.43+00.24 (F) 

18:53:19.1 

01:26:54 

34.43517 

0.24217 

57.1 

3.7 

G2 (MM2) 

G034.77-00.55 (G) 

18:56:50.1 

01:23:11 

34.78117 

-0.56817 

43.6 

2.9 


N otes. Column 1 g i ves th e name of the clump as denoted by Butler & Tan 1 201 2 j) and, in parent hesis, the original name of the clump as assigned 
bv lRathborne et all d2006t ). The F2 clump was not given a designation bv lRathborne et all 120061) . Column 2 gives the name of the IRDC in which 
the clump is embedded. Columns 3 and 4 give the right ascension and declination of the clump and Cols. 5 and 6 give the Galactic longitude and 
latitude of t he clump, as given by lButler & Tanl 1 2012). Column 7 gives the velocity with respect to the local standar d of rest of the parent I RDC, 
as given bv ISimon et all (12006! ). Column 8 gives the kinematic distance o f the cloud from Simon et al. ( 20061). While iKuravama et all 1201 il l find 
a 1.56 kpc distance for IRDC F, containing FI and F2. Foster et al., (2012) find a d istance consistent w ith the kinematic distance for IRDC F, based 
upon extinction measurements. We thus elect to use the kinematic distances from lSimon et al.1 ( 2006 ) for all four clumps in this paper. 

Table 2: Herschel Observation Setup Parameters 


Clump 

Line 

RA 

Dec 

t 

HPBW 

Band 

V 

Eup/k 

ID 

^7mb 

Date 



(h:m:s) 

(°:':") 

(s) 

(arcsec) 


(GHz) 

(K) 



(yy/mm/dd) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 

(12) 

Cl 

8 ->7 

18:42:47.29 

-4:04:06.1 

265.5 

23 

3b 

921.79970 

199 

1342251119 

0.628 

12/09/18 

Cl 

9 -> 8 

18:42:47.29 

-4:04:06.1 

1332 

20 

4a 

1036.91239 

249 

1342254432 

0.635 

12/10/24 

Cl 

ON 

T 

O 

18:42:47.29 

-4:04:06.1 

1425.6 

19 

5a 

1151.98544 

304 

1342253926 

0.590 

12/10/26 

FI 

8 —> 7 

18:53:15.65 

1:26:21.8 

1110 

23 

3b 

921.79970 

199 

1342252136 

0.628 

12/10/03 

FI 

9 —> 8 

18:53:15.65 

1:26:21.8 

1332 

20 

4a 

1036.91239 

249 

1342254434 

0.635 

12/10/25 

FI 

On 

T 

O 

18:53:15.65 

1:26:21.8 

1425.6 

19 

5a 

1151.98544 

304 

1342253937 

0.590 

12/10/26 

F2 

8 —>7 

18:53:19.22 

1:26:47.2 

814 

23 

3b 

921.79970 

199 

1342244947 

0.628 

12/04/25 

F2 

9 —> 8 

18:53:19.22 

1:26:47.2 

1776 

20 

4a 

1036.91239 

249 

1342254433 

0.635 

12/10/25 

F2 

o\ 

T 

O 

18:53:19.22 

1:26:47.2 

1425.6 

19 

5a 

1151.98544 

304 

1342253936 

0.590 

12/10/26 

G2 

8 —> 7 

18:56:49.67 

1:23:15.4 

265.5 

23 

3b 

921.79970 

199 

1342244946 

0.628 

12/04/25 

G2 

9 —> 8 

18:56:49.67 

1:23:15.4 

1332 

20 

4a 

1036.91239 

249 

1342254351 

0.635 

12/11/03 

G2 

10—>9 

18:56:49.67 

1:23:15.4 

1425.6 

19 

5a 

1151.98544 

304 

1342253938 

0.590 

12/11/03 


Notes. Column 1 gives the name of the clump observed while Col. 2 gives the 12 CO transition observed. Columns 3 and 4 are the right ascension 
and declination, respectively, of the map centres. The on source integration time of the observation is given in Col. 5. The average half power 
beam width of the H and V polarizations is given in Col. 6 and the Herschel frequency band of the observation is given in Col. 7. Columns 8 
and 9, respectively, give the res t frequency and upp er level energy, in units of kelvin, of the observed transition from the Spectral Line Atlas of 
Interstellar Molecules (SLAIM:l Remiian et al.ll2007l) . Columns 10 and 11 give the observation ID and the main beam conversion factor used. The 
main beam conversion factor used is the average of the appropriate values for the H and V polarizations based on the revised values of Jellema et 
al. (in prep.). Finally, Col. 12 gives the date the observation was conducted. 


smoothed to match the slightly larger beam size of the CO J 
= 8 —> 7 maps. The resulting data cubes were then converted 
from frequency space into velocity space, with respect to the lo¬ 
cal standard of rest, and then smoothed to a resolution of 2 km 
s 1 (but with 1 km s 1 sampling) from their original 1 MHz reso¬ 
lution (~ 0.3 km s 1 ). To convert between antenna temperatures 
and main beam temperatures, main beam efficiencies of 0.628, 
0.635, and 0.590 were used for the CO / = 8 —* 7, 9 —* 8, and 
10 —♦ 9 data, respectively (Jellema et al., in prep.). A forward 
efficie ncy of 0.96 is used for all observations ( Roelfsema et alJ 
120121) . All intensities quoted in this paper are main beam tem¬ 
peratures, unless otherwise stated. 

To fit a Gaussian profile to each pixel of the data cubes, the 
HIPE multifit graphical user interface (GUI) was used. The pixel 
with the strongest apparent detection was first identified by eye 
and the GUI was used to fit a Gaussian to this spectrum. The 
best fit from this pixel was then used as the initial estimate for 
the automated fitting of the remainder of the pixels. Each fit was 


reviewed and any fits which appeared exceedingly poor were re¬ 
done individually using the GUI. For any spectrum where it was 
ambiguous whether small features near the edges of the line were 
part of the line or just noise, the default fit from the multifit GUI 
was adopted. Since the hot gas component generated from turbu¬ 
lent dissipation should be relatively widespread across a molecu¬ 
lar cloud, a map average spectrum was also created for each data 
set and a Gaussian was fit to these spectra using the HIPE fitting 
GUI. Since the CO J = 10 —* 9 line is not detected in any map, 
the spatially averaged CO J - 10 —* 9 spectra from all four re¬ 
gions were combined in IDL, after first shifting the velocity zero 
points of each spectra. For the Cl, FI, and F2 data, the velocity 
zero point used was chosen to be the average central velocity of 
the spatially averaged CO / = 8 —» 7 and 9 —» 8 spectra. For G2, 
the average of the NiD/_central velocities of the G2-N and G2-S 
cores was used dTan et all2013h . 

While the GUI provides uncertainty estimates for most of 
the fit parameters, we choose to estimate the uncertainty in the 
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integrated intensity from the equation 

dl = RMS x V FWilM x 6v, (1) 

where dl is the uncertainty in the integrated intensity, RMS is 
the root mean square of the baseline, FWHM is the full width 
at half maximum of the line, and 5v is the velocity resolution. 
The RMS of the baseline is determined by applying the HIPE 
statistics task to the regions of the spectrum that appear to be 
devoid of any lines. 

The absolu te flux calibration for these HIFI bands is of the 
order of 10% dRoelfsema et al.ll2012h . but has not been included 
in any below stated intensity uncertainties. The pointing accu- 
racy of Herschel is of the order of 2” dRoelfsema et al.ll2012i) . 

Any spectrum with an integrated intensity less than four 
times the uncertainty in the integrated intensity is classified as 
a non-detection. Spectra with integrated intensities between four 
and five times the uncertainty are classified as tentative detec¬ 
tions and those spectra with integrated intensities greater than 
five times the uncertainty are classified as strong detections. Ten¬ 
tative and strong detections will jointly be referred to as detec¬ 
tions and all further analysis in this paper will consider strong 
and tentative detections equally. We chose a cutoff at four times 
the uncertainty, rather than three times, because too many lines 
with integrated intensities between 3 and 4 times the uncertainty 
in the integrated intensity did not appear to be obvious, good 
quality fits by eye. 

For non-detections, an upper limit to the integrated intensity 
is assigned via 

I xim = 4 RMS x V FWHM me X 6v, (2) 

where FWHM dve is the average full width at half maximum of 
all of the detected lines from the same map. Due to a lack of 
detected lines in the CO J - 10 —> 9 maps, the average FWHM 
of the CO / = 9 —> 8 map of the corresponding region is adopted 
in order to calculate the upper limits of the 10 —» 9 spectra. For 
the G2 CO / = 9 —» 8 and 10 —> 9 maps, where no lines are 
detected, the mean of the FWHM dve values from the Cl, FI, and 
F2 CO / = 9 —> 8 maps (5.4 km s _1 ) is used to estimate the 
integrated intensity upper limit. 

3.2. Results 

Figure [2] shows the spectra and the associated Gaussian fits for 
the three observed lines towards the positions of the Cl-N and 
Cl-S cores. Figure [3] shows the spectra obtained by averaging 
all of the data, for each transition, over the Cl clump and the 
corresponding fits. These spectra are typical of strong detections 
(and non-detections in the case of the 10 —> 9 lines). 

Figures[4]and0show the integrated intensities of the CO J = 
8 —> 7 and 9 —> 8 transitions towards the FI, F2, and Cl clumps. 
Given the close proximity of FI and F2, both clumps are placed 
on the same figure. The G2 clump is not shown, as the CO J = 
8 —» 7 line is only tentatively detected in one location within G2 
and the CO / = 9 —> 8 line is not detected towards any location 
within G2. Maps for the CO J = 10 —» 9 transition are also not 
shown as this line is not detected towards any location. Figure 
[6] shows the integrated intensities of the CO / = 8 —> 7 and 9 
—> 8 transitions towards the FI, F2, and Cl clumps plotted as 
contours overtop of the mass surface densities of these regions, 
for easier comparison between the different data sets. 

The results of the Gaussian fits to the individual pixels are 
given in Table [3] This table lists the fraction of pixels with ten¬ 
tative and strong detections; the maximum peak intensity, inte¬ 


grated intensity, and FWHM; the average peak intensity, inte¬ 
grated intensity, FWHM, and central velocity with respect to the 
local standard of rest; and the minimum FWHM. Table [4] gives 
the best fit parameters for the spatially averaged spectra. Note 
that in Tabled the typical upper limit to the integrated intensi¬ 
ties is given (Ju m , see Eq. ([2} ), which is equal to four times the 
typical integrated intensity uncertainty, while in Tabled it is the 
integrated intensity uncertainty that is instead quoted (dl, see Eq. 
0 ). 

In the Cl, FI, and F2 maps, the CO / = 8 —» 7 line is detected 
in roughly half of the pixels and the 9 —> 8 line is detected in over 
25% of the pixels. In the G2 map, the CO / = 8 —> 7 line is only 
detected in one pixel and in the spatial average. Furthermore, in 
the G2 map, the CO J - 9 —> 8 line is not detected at any location 
and is not detected in the spatial average. The CO / = 10 —> 9 is 
not detected in any of the four clumps, nor in any of the spatially 
averaged spectra of each clump. 

No line is detected in the CO J = 10 —» 9 spectra obtained 
from combining the spatial averages of all four mapped regions, 
with an upper limit on this CO J - 10 —> 9 integrated intensity 
of 0.26 K km s _1 . It should be noted, however, that given the 
differences in the CO / = 9 —> 8 and 8 —> 7 integrated intensities 
between the four regions, in particular between the G2 clump 
and the other three clumps, it is by no means obvious that the 
CO J - 10 —> 9 intensity should be similar between the four 
regions. Thus, the CO / = 10 —> 9 integrated intensity may still 
be slightly larger than 0.26 K km s 1 in any particular clump. 

The Cl, FI, and F2 maps have reasonably similar integrated 
intensity ranges for detections, with the detections of the CO J = 
8 —> 7 and 9 —> 8 lines having integrated intensities ranging from 
0.45 to 2.63 K km s' 1 and from 0.5 to 2 K km s _1 , respectively. 
The average integrated intensity of a detection of the CO J - 8 
—> 7 line varies from 1.11 to 1.31 to 1.05 K km s _1 between the 
Cl, FI, and F2 maps, and the average detection of the 9 —♦ 8 line 
for these three regions changes from 0.67 to 1.14 to 0.75 K km 
s 1 . The integrated intensity of the spatial averages tends to be 
approximately a factor of two lower than the average detection 
from a single pixel and the emission is clearly spatially inhomo¬ 
geneous within each clump. The pixel with both the largest CO J 
— 8 —> 7 and 9 —> 8 integrated intensity lies just to the northeast 
of the 24 micron source in the FI clump. 

Figure |7] shows the ratio of the CO / = 8 —>7 to 9—>8 inte¬ 
grated intensities for all locations with detections in both lines. 
Table 0 gives the minimum, maximum, and average integrated 
intensity ratio for the Cl, FI, and F2 maps. The CO / = 8 —* 7 
to 9 —» 8 integrated intensity ratio ranges from 0.7 to 3.0, with 
a weighted average ratio of 2.0, 1.7, and 1.6 for all the points 
in clumps Cl, FI, and F2, respectively, where both lines are de¬ 
tected. This weighted average is found by weighting the ratios by 
the inverse of their fractional uncertainties, which are typically 
on the order of 20-25%. The ratios of the integrated intensities 
from the spatially averaged WBS spectra are 2.1+ 0.3, 2.6 ± 0.4, 
and 1.9 + 0.3 for clumps Cl, FI, and F2. 

The ratio of the peak intensity of the CO J - 8 —> 7 and 9 —* 8 
lines are reasonably similar, with weighted averages of 2.1, 2.0, 
and 1.4 for all of the points with detections in both transitions 
and ratios of 2.1 ± 0.5, 2.4 ± 0.5, and 2.4 ± 0.6 for the spatially 
averaged WBS spectra (for Cl, FI, and F2 respectively). 

In the FI, F2, and G2 maps, detection of the CO J = 9 —> 
8 line is always accompanied with a CO / = 8 —> 7 detection; 
there is one location (18 /! :53"T8.5\ 1°:27’:10”) in the F2 map 
where the CO J - 9 —> 8 integrated intensity is larger than the 
CO / = 8 —♦ 7 integrated intensity, due to the 9 —> 8 line having 
a larger FWHM. At this location, the CO J = 8 —> 7 line still has 
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Table 3: Fitting Results for the Herschel Data 


Clump 

(i) 

Line 

(2) 

/tent 

(3) 

./strong 

(4) 

GlB .max 

(K) 

(5) 

T MFS.ave 

(K) 

(6) 

/max 

(K km s~ 
(7) 

lave 

*) (K km s~ 
(8) 

him 

') (KkmC 
(9) 

FWHM max 
) (km s^ 1 ) 
(10) 

FWHM mm 
(km s _1 ) 
(11) 

FWHM me 
(km s _1 ) 
(12) 

Vlsr 
(km s^ 1 ) 
(13) 

Cl 

8 —> 7 

4(10.0%) 

13 (32.5%) 

0.58 

0.27 

2.03 

1.11 

0.67 

10.2 

2.1 

4.6 

79.2 

Cl 

9 —> 8 

7(17.5%) 

8 (20.0%) 

0.21 

0.14 

0.96 

0.67 

0.50 

5.9 

3.5 

4.6 

79.4 

Cl 

10 —* 9 

0 

0 





<1.75 





FI 

8 —> 7 

5 (12.5%) 

19 (47.5%) 

0.41 

0.23 

2.63 

1.31 

0.62 

8.7 

2.3 

5.4 

56.3 

FI 

9 —> 8 

3 (7.5%) 

7(17.5%) 

0.21 

0.15 

2.03 

1.14 

0.75 

9.5 

5.1 

6.9 

55.5 

FI 

10 —* 9 

0 

0 





<2.16 





F2 

8 —* 7 

5 (12.5%) 

18(45.0%) 

0.33 

0.19 

1.96 

1.05 

0.67 

7.4 

3.0 

5.4 

58.0 

F2 

9 —> 8 

3 (7.5%) 

8 (20.0%) 

0.26 

0.15 

1.16 

0.75 

0.54 

5.8 

3.8 

4.8 

57.7 

F2 

10 —* 9 

0 

0 





<2.12 





G2 

8 —* 7 

1 (2.9%) 

0 (0%) 

0.18 

0.18 

0.72 

0.72 

0.64 

3.8 

3.8 

3.8 

62.0 

G2 

9 —> 8 

0 

0 





<0.60 





G2 

10 —* 9 

0 

0 





<1.95 






Notes. Column 1 gives the name of the clump observed while Col. 2 gives the CO transition observed. Columns 3 and 4 give the number of 
pixels in which there are tentative and strong detections, respectively. The number in parentheses is the corresponding percentage of pixels in the 
entire map. The maximum peak intensity and mean peak intensity of all of the detections are given in Cols. 5 and 6, while the maximum and 
mean integrated intensities of all of the detections are given in Cols. 7 and 8, respectively. Column 9 gives four times the average uncertainty 
in the integrated intensities of all of the detections, which gives a measure of the typical upper limits for any non-detections. For maps with no 
detections. Col. 9 gives four times the average uncertainty of all of the pixels. Columns 10-12 give the maximum, minimum, and mean FWHM of 
all detections, respectively. Column 13 gives the mean central velocity with respect to the local standard of rest of all of the detections. 


Table 4: Spatially Averaged Fit Parameters 


Clump 

(1) 

Line 

(2) 

Tmb 

(K) 

(3) 

I 

(K km s' 1 ) 
(4) 

FWHM 
(km s' 1 ) 
(5) 

Flsr 
( km s' 1 ) 
(6) 

Cl 

8 —> 7 

0.14(0.02) 

0.71 (0.05) 

4.8 (0.6) 

79.5 (0.2) 

Cl 

9 —> 8 

0.065 (0.012) 

0.34 (0.04) 

4.9 (1.0) 

78.9 (0.4) 

Cl 

10 —> 9 


...(0.12) 



FI 

8 —> 7 

0.13(0.02) 

0.92 (0.08) 

6.7 (0.9) 

55.6 (0.4) 

FI 

9 -> 8 

0.054 (0.009) 

0.35 (0.04) 

6.1 (1.4) 

55.9 (0.6) 

FI 

Os 

T 

O 


...(0.15) 



F2 

8 —> 7 

0.14(0.02) 

0.69 (0.08) 

4.7 (0.9) 

57.9 (0.4) 

F2 

9 —> 8 

0.056 (0.010) 

0.36 (0.04) 

6.0 (1.3) 

58.4 (0.5) 

F2 

10 —> 9 


...(0.12) 



G2 

8 —>7 

0.05 (0.01) 

0.31 (0.06) 

5.7 (1.6) 

59.7 (0.7) 

G2 

9 8 


... (0.04) 



G2 

Os 

T 

O 


...(0.14) 




Notes. Column 1 gives the name of the clump observed while Col. 2 gives the CO transition observed. Columns 3-6 give the best fit parameters 
for the spatially averaged spectrum from each map. The parentheses give the uncertainties for each parameter, as given by the Gaussian fitting 
routine, except for the integrated intensity, where the uncertainty is given by Eq. CD. Column 3 gives the peak intensity while column 4 gives the 
integrated intensity. The FWF1M is given in Col. 5. The central velocity with respect to the local standard of rest is given in Col. 6. 


a larger peak intensity. There are four locations in the C1 map 
where the CO / = 9 —* 8 line is detected but the CO / = 8 —* 1 
line is not, despite reasonably similar RMS values between the 
two lines. For three of these locations, the CO / = 9 —> 8 line is 
only tentatively detected. The exception is the pixel centered at 
(18 / ':42 m :47.3 s , -4°:04’:21”), where the CO J = 9 -> 8 line has a 
strong detection. Given the generally low signal to noise of all of 
the detected lines, it is possible that these unusual line ratios are 
just due to noise, although the CO J — 9 —> 8 line may indeed 
be stronger than the CO J — 8 —> 7 line towards some, or all, of 
these locations, in particular the Cl location with a strong 9 —> 8 
detection. 

The average linewidth of detections ranges between 4 and 7 
km s' 1 , depending upon the line and the clump. There are no ob¬ 
vious, significant trends in the distribution of FWHM or central 
velocity in the four clumps, although the determination of these 


line properties is significantly limited by the low signal to noise 
of the data. Maps of these quantities are therefore not shown. 
No lines other than the CO lines are detected in these Herschel 
observations. 

In the Cl, FI, and F2 maps, the central velocities of the de¬ 
tected lines are all within 3.5 km s' 1 of the expected velocity, 
with the most discrepant velocities typically being associated 
with either weaker detections or broader lines. For the one pixel 
from the G2 map where the CO / = 8 —> 7 line is detected and 
the spatially averaged CO / = 8 —♦ 7 spectrum of G2, however, 
the central velocities of the detected lines vary from the known 
velocity of the G2 clump by over 15 km s -1 . As such, we con¬ 
sider the G2 detections to be false positives, possibly caused by 
random noise in the spectra. The velocity discrepancy between 
the line detected in the G2 spatial average and the known G2 
central velocity is unlikely to be due to an outflow, as protostel- 
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lar outflows are typically well collimated, such that the outflow 
emission would only contribute to a few pointings within the sur¬ 
veyed area. In such a case, a detection in the G2 spatial average 
would require detections in multiple pixels, since map positions 
off of the outflow axis would contribute no signal to the spatial 
average. We also consider it highly unlikely that this detection 
is a line from a more complex molecular species, given that CO 
should be by far the most abundant molecule in these IRDCs. If 
the G2 detections are false positives due to noise, it is, however, 
unusual that the G2 CO / = 9 —> 8 spatially averaged WBS and 
HRS spectra show an emission feature at a level of three times 
the integrated intensity uncertainty, just below our threshold for 
a detection, at a similar 60 km s -1 central velocity to the CO J = 
8 —> 7 tentative detections, as the central velocities of noise in¬ 
duced false positives in different data sets should not necessarily 
be correlated. 

4. Comparison PDR model 

4.1. Setup 

The [Kaufman et al.l (119991 ) PDR model calculates the emission 
from a semi-infinite, uniform density slab of material that is il¬ 
luminated from one side. The PDR model includes molecular 
freeze out and starts with roughly Solar chemical abundances. 

A microturbulent Doppler broadening of [5 — 1.5 km s [ is 
used, which corresponds to a FWHM of approximately 2.5 km 
s 1 , since FWHM = 2/1Vln(2). This is smaller than the FWHM 
observed for the mid-/ lines, but consistent with the linewidths 
of the previously observed l2 CO / = 3 —> 2 line towards these 
clumps dSanhueza et al.tf2010t Pon et al. in prep., Paper II). Sin¬ 
gle dish observations of the N 2 D + / = 2 —» 1 line, which have 
a similar beam size to that of our Herschel observations, also 
reveal linewidths of the order of 2 km s _I dFontani et al.ll201 ll) . 
Using a smaller FWHM for the PDR model will increase the 
emission in the mid-/ lines, as the lower lines will be less effec¬ 
tive coolants. 

To give the PDR models the best chance at matching the 
bright emission detected in the mid-/ CO fines, we use th e larger 
PAH abundance introduced by iHollenbach et al.l (12012|) . as this 
enhanced PAH abundance produces extra heating in the periph¬ 
ery of a cloud which leads to much larger mid-/ CO fluxes. 

Th e six cores surveyed have densities of the order of 10 5 
cm' 3 dButler & Tanl l2009i l20ll iTan et all l2()H iButler et all 
[20f4l) . but the gas surrounding these cores is likely at a 

Table 5: CO / = 8 —>7 to 9—>8 Integrated Intensity Ratios 


lower density closer to 10 4 cm 3 , as t ypically found in IRDCs 
(iRathborne et alj2006fclTan et al.l2014l) . Most of the higher / CO 
emission is also expected to come from the warmer, outer layers 
of a molecular cloud, where the density is lower than in the cen¬ 
ters of dense cores. We thus compare our Herschel intensities to 
PDR models with densities of either 10 4 or 10 5 cm -3 . 

We examine models with interstellar radiation fields (ISRFs) 
of log(Go) = 0 (Go = 1 Habing), where the average far-ultraviolet 
ISRF in free space is 1.7 Habing, or 1.6 x ICC 3 erg ctrT 2 s _1 . We 
also examine a model with a density of 10 4 cm -3 and an ISRF 
of log(Go) = 0.5 (Go » 3 Habing), to test the effect of increasing 
the UV field. 

In all of the PDR models, the CO / = 8 —* 7 and higher lines 
are optically thin (r < 0.5), such that we double the integrated 
intensities of the lines to account for emission directed inwards 
from the far side of the cloud. Since the models go to an Ay of 
10.4 mag, this doubling accounts for the emission from a slab 
with an Ay of 20.8 mag. While the visual extinctions towards 
the IRDC core centers are much higher than this, of the order 
of 100 mag, we consider the PDR model predictions to still be 
reasonably accurate, as most of the mid-/ CO emission comes 
from the warm outer layers where CO is not heavily frozen out. 
In the three models, over 90% of the emission in the CO / = 8 
—> 7, 9 —» 8, and 10 —> 9 lines is emitted from gas at Ay < 1.5. 

4.2. PDR results 

Figure [8] shows the PDR model predictions as well as the ob¬ 
served integrated intensities of the mid-/ CO lines. All three 
PDR models underpredict the observed CO / = 9 —> 8 integrated 
intensities in the Cl, FI, and F2 clumps, including the integrated 
intensities from the spatial averages, which have almost an or¬ 
der of magnitude lower integrated intensity than the maximum 
CO / = 9 —> 8 integrated intensity measured from an individual 
pixel. The models with an ISRF of 1 Habing underpredict the 
observed CO / = 8 —> 7 integrated intensity, although the higher 
ISRF model for a density of 10 4 cm 3 is capable of producing an 
8 —» 7 integrated intensity larger than the largest observed inte¬ 
grated intensity of the 8 —> 7 line. The non-detection of the CO 
/ = 10 —> 9 line in all of the spectra is consistent with all three 
PDR models, while the non-detection of any mid-/ CO lines in 
the G2 clump is only consistent with the highest density PDR 
model. The 0.26 K km s _I upper limit on the CO / = 10 —> 9 in¬ 
tegrated intensity from the combined spectrum of all four clumps 
is still consistent with all three PDR models. 


Clump 

(1) 

/ratio 

(2) 

min 

(3) 

max 

(4) 

ave 

(5) 

WBS 

(6) 

Cl 

11 (27.5%) 

1.2 

2.9 

2.0 

2.1 (0.3) 

FI 

10 (25.0%) 

1.1 

2.6 

1.7 

2.6 (0.4) 

F2 

11 (27.5%) 

0.7 

3.0 

1.6 

1.9 (0.3) 


Notes. Column 1 gives the name of the clump observed. Column 2 gives 
the number of pixels where both the CO / = 8 —» 7 and 9 —» 8 lines are 
detected while the value in parentheses gives the corresponding fraction 
of pixels with detections in both lines. Columns 3-5 give, respectively, 
the minimum, maximum, and weighted average ratio of the CO / = 8 —> 
7 to 9 —» 8 integrated intensities. The integrated intensity ratio of these 
two lines from the spatially averaged spectra is given in Col. 6, with the 
number in parentheses being the uncertainty in the ratio. 
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The log(n) = 4, log(Go) = 0 model produces a large CO / = 8 
—♦ 7 to 9 —♦ 8 integrated intensity ratio of 13.6, while the higher 
ISRF model for this density produces a slightly lower ratio of 
9.2. The log(n) = 5, log(Go) = 0 model generates a ratio of 7.4. 
For all three models, the ratio of the CO / = 8—>7to9—>8 
integrated intensities is nowhere near the average observed ratios 
of 1.6 to 2.0 for individual pixels in the Cl, FI, and F2 maps. 
The models are also highly inconsistent with the CO / = 8 —> 7 
to 9 —♦ 8 integrated intensity ratios obtained from the spatially 
averaged spectra, which are in the range of 1.9 to 2.6 for the Cl, 
FI, and F2 maps. 

5. DISCUSSION 

Based on the PDR models, the CO / = 9 —* 8 line should not 
have been detected in any pixel or in any of the spatially aver¬ 
aged spectra. Instead, the observations clearly show CO / = 9 
—> 8 detections in over 25% of the pixels in the Cl, FI, and F2 
maps, as well as in the spatially averaged spectra of these maps. 
Similarly, the PDR models predict a much larger CO / = 8 —> 7 
to 9 —♦ 8 integrated intensity ratio than observed, with the dis¬ 
crepancy between the highest observed ratio calculated from the 
spatially averaged spectra and the lowest predicted ratio still be¬ 
ing 12 times the uncertainty for the observed ratio. While the 
average ratio from pixels with detections in both lines may be 
slightly skewed to lower ratios, since requiring both lines to be 
detected preferentially selects pixels with larger CO / = 9 —> 
8 integrated intensities, the spatially averaged spectra, which do 
not have this selection bias, produce similarly low ratios for the 
Cl, FI, and F2 clumps. 

Increasing the density of the PDR does not significantly 
change the absolute integrated intensities of the mid-/ lines, but 
rather primarily decreases the ratio of the CO / = 8 —>7 to 9—> 

8 integrated intensities. Raising the ISRF increases the intensi¬ 
ties of all transitions and slightly lowers the CO / = 8 —> 7 to 

9 —> 8 ratio. For relatively weak ISRFs (Go < 100), increasing 
the ISRF has a very limited effectiveness in creating hot CO, as 
larger ISRFs increase the size of the bounding C + layer, causing 
CO to form at larger extinctions, where the heating rate, and thus 
temperature, are lower. CO is also such an effective coolant, that 
upon formation of significant quantities of CO, the gas tempera¬ 
ture tends to drop rapidly. This suggests that a PDR region with 
a much higher density and ISRF would be needed to consistently 
match both the absolute integrated intensities of the CO / = 8 —» 

7 and 9 —» 8 lines and the ratio between these lines. Given the 
observational constraints on the density and ISRF of these re¬ 
gions (see Sect. 0. it seems highly unlikely that PDR emission 
alone can explain the observed mid-/ CO line intensities. The 
inability of the PDR models to reproduce the large observed CO 
/ = 9 —> 8 integrated intensities and the CO / = 8 —>7 to 9 —* 

8 integrated intensity ratios strongly suggests that there is a gas 

component in the Cl, FI, and F2 clumps that is not accounted 
for by the PD R models. __ 

RADEX (Ivan der Tak et al.l 12007 ) ID, non-local thermody¬ 
namic equilibrium (LTE) radiative transfer models with CO col¬ 
umn densities of the order of 2 x 10 18 cm -2 , line widths of 5 
km s -1 , densities of 10 5 cm -3 , and gas temperatures around 20 
K also produce CO / = 8 —>7 to 9—>8 intensity ratios of the 
order of 10. Only at temperatures above 40 K does the line ratio 
approach the observed values, although such high temperature 
models, with the full CO column density of 2 x 10 18 cm' 2 , sig¬ 
nificantly overproduce, by an order of magnitude, the observed 
integrated intensities. Radex models with lower column densi¬ 
ties, such that the CO / = 8 —» 7 and 9 —» 8 lines are optically 


thin, suggest a ratio between the lines of approximately 2 is only 
reached at temperatures above 100 K. Thus, the low observed ra¬ 
tio of the CO / = 8 —>7 to 9—>8 lines suggests that the missing 
gas component in the PDR models is a hot gas component, with 
a temperature of at least 40 K, and likely above 100 K. 

The variation of the mid-/ CO line strengths and ratios 
across the surveyed regions implies that the generating source 
of the hot gas component must be spatially intermittent and can¬ 
not be a relatively uniform heating source, such as cosmic ray 
heating. The beam size changes by less than 10% between the 
CO / = 8 —> 7 and 9 —> 8 transitions, such that differing beam 
filling factors due to changes in the beam size cannot account 
for the line ratio differences between the observations and the 
models. 

The dissipation of turbulent motions in low velocity shocks 
is expected to produce a hot gas component filling a small vol¬ 
ume fraction of the cloud, and such a hot gas component has 
already been observed towards a low mass star-forming region, 
Perseus B1-E5 ( Pon et all 12012i 12014i) . The large mid-/ CO in¬ 
tegrated intensities towards the Cl, FI, and F2 regions may be 
due to such a low volume filling factor hot gas component gen¬ 
erated by the dissipation of turbulence. Preliminary models of 
turbulence dissipating in low velocity shocks propagating in 10 5 
cm' 3 gas predict CO / = 8 —» 7 and 9 —> 8 integrated intensi¬ 
ties on the order of 1.5 and 1.0 K km s _1 , respectively, which 
are comparable to the observed intensities in the Cl, FI, and F2 
clumps (Pon et al., in prep.. Paper III). 

While the Cl-N, Cl-S, FI, and F2 cores are believed to be 
starless, there are more evolved sources within IRDCs C and 
F that could be responsible for generating a warm gas compo¬ 
nent. For instance, as shown in the middle panel of Fig. Q] there 
is a 24 micron source towards the northeast corne r of the FI 
clump (Sheph erd et alj|2007t IChambers et akll2009l) . indicative 
of an embedded young stellar object. iFoster et al. (1201 4li have 
also identified an embedded population of low mass protostars 
near this 24 micron source that could heat gas from protostellar 
outflows. This 24 micron source location is roughly coincident 
with the largest CO / = 8 —> 7 and 9 —» 8 integrated intensities 
detected in any map. The mid-/ CO emission in the F2 map also 
preferentially comes from the vicinity of the MM7 clump, which 
is known to contain a water maser , a s ig npo st of embedded high 
mass star formation (IChambers et al.ll2009t) . It is less obvious, 
however, what embedded sources, if any, could be responsible 
for the mid-/ emission in the Cl map. The average FWHM of 
the mid-/ CO lines is of the order of 5 km s _1 , which is slightly 
larger than the 3 km s _I line width of the 12 CO / = 3 — * 2 line 
seen towards these clumps (ISanhueza et al.ll2010t Paper II). The 
larger linewidths of the higher lines could be indicative of broad¬ 
ening from outflows, but could also be due to the low signal to 
noise of the mid-/ detections. High velocity outflows generate 
very strong shocks that can lead to gas temperatures above 1000 
kelvin (e.g JKaufman & NeufelcM1996h . which would create dif¬ 
ferent line ratios compared to lower velocity shocks. 

While the ratios of the mid-/ CO lines can discriminate be¬ 
tween different origins for the hot gas component by constrain¬ 
ing the temperature of the hot gas. Fig. [8] shows that the CO / 
= 8 —» 7 emission could come from either the hot gas compo¬ 
nent or the ambient gas, depending upon the density of the gas 
and the ISRF. Observations of lower / CO lines are thus vital for 
constraining the PDR parameters and determining what fraction 
of the 8 —> 7 emission comes from the hot gas component. Ob¬ 
servations of lower / CO transitions towards the Cl, FI, and F2 
clumps will be presented in a future paper (Paper II). The sim¬ 
ilarity in the spatial distribution of the CO / = 8 —♦ 7 and 9 —» 
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8 emission, however, suggests that the two lines are powered by 
the same excitation process, such that the 8 —> 7 emission may 
also come primarily from the warm gas component. 

The G2 region is clearly different from the Cl, FI, and F2 
regions, in that no emission in any of the observed lines is de¬ 
tected in the G2 maps. It is not obvious, however, whether there 
are significant differences between the Cl, FI, and F2 regions, 
given the large variance in emission between pixels in each of 
the individual maps. The lack of detections in the individual pix¬ 
els in G2 can be attributed to the lack of a hot gas component 
within G2, or at least a less substantial hot gas component in 
G2 than in the other three observed regions. Of the four clumps 
observed, the G2 clump has the smallest measured N 2 FC line 
widths dFontani et al.l 120171) . such that G2 appears to be less 
turbulent than the other clumps. This lower level of turbulence 
may account for a reduced hot gas component, if the hot gas 
components in the Cl, FI, and F2 clumps are due to turbulent 
dissipation. The lower turbulent state of G2 may be related to 
its evolved state, given that a transition to coherence is seen in 
move evolved, lower mass prestellar cores (iPineda et al.ll20ldl) . 
The G2 clump also has almost twice the D/H ratio of the other 
three clumps as measured from the ratio of N 2 D + and N 2 FC 
dFontani et al.ll201 Tl) and an orde r of magnitude large r D/H ratio 
measured from NH 3 and NH 2 D dFontani et al.ll2014h . thus sug¬ 
gesting that the G2 clump is cooler and potentially more evolved. 
Alternatively, if G2 is younger, then G2 would be more likely 
to be free of embedded low mass protostars and the associated 
heating from these protostars and their protostellar outflows. 

While the non-detection of the CO / = 8 —> 7 line in the spa¬ 
tial average of the spectra from the G2 map is consistent with the 
highest density PDR model, the G2 clump is believed to have the 
lowes t density of the four clumps observed (But ler & Tanl l2009. 
1 201 2l) . A lower PAH abundance than assumed in the PDR mod¬ 
els would lead to cooler temperatures at low visual extinctions, 
which would in turn reduce the expected integrated intensities 
of the mid-/ CO lines from the PDR zone. As such, a lower 
PAH abundance in the G2 clump would allow the G2 clump non¬ 
detection to be still consistent with the expected CO / = 8 —> 7 
emission from a lower density PDR zone around the G2 clump. 

6. CONCLUSIONS 

Maps of the CO / = 8 — * 7, 9 —» 8 , and 10 — * 9 transitions 
were obtained with the Herschel Space Observatory towards the 
Cl, FI, F2, and G2 clumps. These four clumps are all embedded 
within IRDCs and each contains at least one core that appears to 
be quiescent and prestellar. The CO / = 8 —> 7 and 9 —> 8 lines 
are detected in roughly 25 to 50% of the pixels towards Cl, FI, 
and F2, but are not reliably detected in any pixel or in the spatial 
average of all pixels towards the G2 clump. The CO / = 10 —» 

9 line is not detected in any pixel or any spatial average, thereby 
placing a 0.5 K km s 1 upper limit on the integrated intensity in 
the spatially averaged spectra. The detected CO J - 8 —> 7 and 
9 —» 8 lines have typical integrated intensities of the order of 1.2 
and 0.9 K km s ' 1 in individual pixels and integrated intensities 
of the order of 0.8 and 0.35 K km s 1 in the spatially averaged 
spectra. The emission is, however, not spatially uniform, with 
the line emission varying by at least a factor of four between 
different regions within the same map. The average integrated 
intensity ratio of the CO / = 8 —>7 to 9 —>8 lines is found to be 
of the order of 1.6 to 2.6 

PDR models based on lKaufman et al.l (119991 ) were compared 
to the observed integrated intensities. While the CO / = 8 —> 7 
integrated intensities can be explained by PDR models with rea¬ 


sonable densities and ISRFs, the models all underpredict the CO 
J - 9 —> 8 integrated intensity, with the models predicting that 
the CO J - 9 —> 8 line should have been undetected at all loca¬ 
tions given the RMS level obtained. The PDR models also pro¬ 
duce CO / = 8 —>7 to 9 —>8 integrated intensity ratios between 
7 and 14, which are significantly larger than the observed ratios. 
The combination of abnormally high integrated intensities of the 
CO J - 9 —> 8 line and low ratios between the CO / = 8 —> 7 and 
9 —> 8 lines strongly suggests the presence of a hot gas compo¬ 
nent within the Cl, FI, and F2 clumps with a temperature likely 
above 100 K. The presence of a hot gas component is consistent 
with predictions for the dissipation of turbulence in low velocity 
shocks dPon et al.ll2012 ). but could also be due to heating from 
outflow generated shocks or embedded protostars. The emission 
around the FI and F2 clumps, in particular, is concentrated near 
regions with known star formation, although there are no obvi¬ 
ous embedded sources associated with the emission regions near 
the Cl clump. 
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Fig. 1: Mass surface density derived bv i Butler & Tanl ( 2012 '). The top, bottom left, and bottom right panels show IRDCs F, C, and 
G, respectively. The contours start at 0.075 g cm -2 (Ay of 17 mag) and increase by increments of 0.075 g cm~ 2 . The blue diamonds 
give the central locations of the cores as seen in NVD 1 emission dTan et alil2013l) . The regions shown in the bottom two panels 
(IRDCs C and G) are the regions mapped by Herschel in the CO / = 9 —> 8 line. The large blue rectangles in the top panel (IRDC 
F) show the areas mapped by Herschel in the CO / = 9 —> 8 line for the FI and F2 clumps. In the bottom left panel , the mass surface 
density hole near 18 , ':42"':48.0 V , -4°:04’:00” is due to the presence of an 8 micron source at this location preventing^ reliable mass 
surface density from being determined from the extinction of the 8 micron Galactic background (IButler et all 120141) . Such mass 
surface density holes, corresponding to locations of other 8 micron sources, can be seen in the other maps. In the top panel , the 
cross gives the location of the 24 micron source associated with the FI clump ((Chambers et al.ll2009l) and the “X” gives the location 
of the lRathborne et all (120061) MM7 clump. 
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Fig. 2: CO / = 8 —* 7, 9 —> 8, and 10 —* 9 lines (left, middle, right) towards the Cl-N (top) and Cl-S ( bottom ) cores. The best fitting 
Gaussians are overlaid in red. The integrated intensity of the fit and its uncertainty are given in K km s -1 in the top right corner of 
each spectral box. For non-detections, the upper limit in K km s _1 is instead given. 





Velocity (km s - ') Velocity (km s _l ) Velocity (km s - ') 

Fig. 3: Spatially averaged spectra of the CO / = 8 —* 7, 9 —* 8, and 10 —> 9 lines (left, middle, right) towards the Cl clump. The 
best fitting Gaussians are overlaid in red. The integrated intensity of the fit and its uncertainty are given in K km s 1 in the top right 
corner of each spectral box. For non-detections, the upper limit in K km s 1 is instead given. 
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Fig. 4: Integrated intensities of the CO / = 8 —» 7 (l eft) and CO / = 9 —* 8 (right) lines towards the Cl clump are shown in the colour 
scale. The contours are mass surface density from lButler & Tanl d2012l) . with the contours starting at 0.075 g cm~ 2 and increasing 
by increments of 0.075 g crcT 2 . The blue diamonds are the locations of the Cl-N and Cl-S cores, with the top left core being Cl-N. 
The blue circles have diameters equal to the half power beam width of the Herschel beams. 
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Fig. 5: Integrated intensities of the CO / = 8 —» 7 (left) and CO /_= 9 —* 8 (right) lines towards the FI and F2 clumps are shown 
in the colour scale. The contours are mass surface density from iButler & Tatil d2012l) . with the contours starting at 0.075 g cm~ 2 
and increasing by increments of 0.075 g cm -2 . The blue diamonds are the locations of the FI and F2 cores, with the top left core 
being F2. The large blue rectangles denote the area surveyed with Herschel. The cross gives the location of the 24 micron source 
associated with the FI clump (IChambers et al.120091 ) and the “X” gives the location of the Rathborne et aD (i2006t) MM7 clump. The 
blue circles have diameters equal to the half power beam width of the Herschel beams. 



Right Ascension (J2000) Right Ascension (J2000) 

Fig. 6: Integrated intensities of the CO J = 8 —> 7 (green) and CO J = 9 —» 8 (purple) lines are shown as the thick contours. The left 
panel shows the Cl clump while the right panel shows the FI and F2 clumps. The contours start at four times the average uncertainty 
and increase in increments of two times the average uncertainty (see Col. 9 of Table [3}. For the right panel, the average uncertainty 
for each line was taken as the average value of between the FI and F2 maps. The colour scale and thin white contours are mass 
surface density from lButler & Tant ( 2012 ). with the contours starting at 0.075 g cm~ 2 and increasing by increments of 0.075 g cm~ 2 . 
The blue diamonds are the locations of the Cl -N, Cl-S, FI, and F2 cores. The cross gives the location of the 24 micron source 
associated with the FI clump (IChambers et al.ll2009l) and the “X” gives the location of the iRathborne et ah! d2006l) MM7 clump. 
The large blue rectangles on the right panel denote the area surveyed with Herschel, while the left panel shows the entire Herschel 
survey region around the Cl clump. 
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Fig. 7: Ratios of the integrated intensities of the CO / = 8 —>7 to 9—>8 lines for locations where both lines are detected are shown 
in the colour scale. The left panel shows the Cl clump while the right panel shows the FI and F2 clumps. The contours are mass 
surface density from' Butler & Tani d2012l) . with the contours starting at 0.075 g cm~ 2 and increasing by increments of 0.075 g cm~ 2 . 
The blue diamonds are the locations of the Cl-N, Cl-S, FI, and F2 cores. See Fig.[T|for the labels of these cores. The cross gives 
the l ocati on of the 24 micron source associated with the FI clump (IChambers et alJl2009t> and the “X” gives the location of the 
iRathborne et alJ (120061) MM7 clump. The large blue rectangles denote the area surveyed with Herschel. 
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Fig. 8: Integrated intensities of various CO rotational transitions as predicted by the lKaufman et al. ( 1999 ) PDR models (lines). The 
solid line is for a PDR model with log(n) = 4 and log(Go) = 0, the dotted line is for log(n) = 4 and log(Go) = 0.5, and the dashed line 
is for log(u) = 5 and log(Go) = 0. The black diamonds give the integrated intensities from individual pixels. The light blue points 
are centered on the mean value of all of the detections towards a particular clump and the error bars on these points are sized to 
the average uncertainty of all single pixel detections, such that the uncertainty in the average integrated intensity detection is much 
smaller. For maps with no detections, the light blue upper limit instead shows the average upper limit for an individual pixel. The 
red points show the integrated intensity of the average of all of the spectra across each map, with the error bars showing four times 
the uncertainty of this integrated intensity. The top left, top right, bottom left , and bottom right panels show the data for the Cl, FI, 
F2, and G2 regions, respectively. Note how none of the PDR models can account for the large observed CO / = 9 —> 8 integrated 
intensities. 
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